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Addressing the question of whetlier the Hawking effect depends on degrees of freedom at ultra- 
high (e.g., Planckian) energies/momenta, we propose three rather general conditions on these degrees 
of freedom under which the Ifawking effect is reproduced to lowest order. As a generalization of 
Corley's results, we present a rather general model based on non-linear dispersion relations satisfying 
these conditions together with a derivation of the Hawking effect for that model. However, we also 
demonstrate counter-examples, which do not appear to be unphysical or artificial, displaying strong 
deviations from Hawking's result. Therefore, whether real black holes emit Hawking radiation 
remains an open question and could give non-trivial information about Planckian physics. 
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I. INTRODUCTION 

The striking similarity between the laws of black hole 
physics and the (zeroth till third) law of thermodynam- 
ics motivated the idea to assign thermodynamic prop- 
erties such as temperature and entropy to black holes 
0. Hawking's prediction 0] that black holes should 
emit thermal radiation with the temperature being con- 
sistent with the thermodynamic interpretation strongly 
supported this idea. As a consequence, the concept of 
black hole entropy as given by the surface (horizon) area 
of the black hole in Planckian units (instead of the vol- 
ume, for example) is now used in many ways to estimate 
the total entropy of other objects - which is expected 
to be a measure of the number of fundamental degrees 
of freedom of the underlying theory (including quantum 
gravity). 

However, in view of the (exponential) gravitational 
red-shift near the horizon, the outgoing particles of the 
Hawking radiation originate from modes with extremely 
large (e.g., trans-Planckian) wavenumbers. As the known 
equations of quantum fields in curved space-times are ex- 
pected to break down at such wavenumbers, the deriva- 
tion of the Hawking radiation has the flaw that it applies 
a theory beyond its region of validity. This observation 
poses the question of whether the Hawking effect is inde- 
pendent of Planckian physics or not. 

One way to address this question is to model the break- 
down of the (usual) local Lorentz invariance (to be ex- 
pected at the Planck scale) by a (non-linear) deviation 
from the linear dispersion relation at high wavenumbers, 
see, e.g., This method is inspired by the black hole 

analogues which exploit the analogy between the propa- 
gation of excitations (e.g., sound waves) in laboratory- 
physics systems and quantum fields in curved space- 
times, see, e.g., HIEH- 

In Sections m IIIII IIVI and we generalize and sim- 
plify the model and the results presented by Corley in 3] 
(see also 0) trying to identify and to present the crucial 
points. Section IVII is devoted to the question of which 



conditions and assumptions regarding Planckian physics 
are needed to reproduce Hawking's result - together with 
some counter-examples. 

II. LINEAR MODEL 

At first we consider a sub-luminal dispersion relation, 
cf. Fig. ^ which is in some sense conceptually more clear 
because the in-modes generating the Hawking radiation 
come from outside the black hole. The horizon acts as 
a classical turning point where the JWKB (geometric 
optics) approximation breaks down allowing phenomena 
like particle creation. In contrast to Ref. 3] , we shall not 
specify the shape of the dispersion relation apart from 
some rather general assumptions. 

A. Wave Equation 

The geometry as seen by the low-energy particles is 
described in terms of the 1-1-1 dimensional Painleve- 
GuUstrand-Lemaitre @ metric {h = c = 1 throughout) 

ds^ = - [dx - v{x) dtf 

= [I- v'^]dt^ + 2vdtdx ~ dx^ . (1) 

The quantity v{x) can be interpreted as the local veloc- 
ity of the freely falling frames measured with respect to 
the time t corresponding to the Killing vector dt of that 
stationary metric. In terms of the sonic black hole ana- 
logues, t is the laboratory time and v is just the position- 
dependent velocity of the fluid with the (assumed to 
be constant) speed of sound being absorbed by a re- 
definition of the coordinates. Since the behavior near 
the horizon in arbitrary dimensions is essentially 1-1-1 di- 
mensional for each mode, we restrict ourselves to 1+1 
dimensions. Furthermore, we neglect back-scattering (as 
induced by the angular-momentum barrier, for example). 

In order to ensure hyperbolicity, causality, and stabil- 
ity, we only allow second time-derivatives. Hence the 
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generalized Klein-Fock-Gordon equation reads 
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with the function F representing the non-trivial disper- 
sion relation. In general, the function F might contain an 
arbitrary number of derivatives - i.e., be non-local (think 
of a lattice, for example). Note that we do not take into 
account absorption (i.e., F is real). The resulting dis- 
persion relation (w -I- vk)"^ = fc^ — F[— /c^] is plotted in 
Fig. [2 
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FIG. 1: Sub-luminal dispersion relation (not to scale). The 
points of intersection (black circles) with the two lines for 
w > 1 (i.e., a; < 0) and v < 1 (i.e., i > 0) determine the solu- 
tions of the dispersion relation for a given oj. The two points 
corresponding to large wavenumbers have group velocities 
smaller than 11, i.e., they are "swept away" and approach the 
horizon from above a; > 0. Hence these solutions are the in- 
modes. The other solution at a; > with the group velocity 
exceeding v represents the outgoing Hawking radiation. The 
only solution beyond the horizon a; < again has a group ve- 
locity smaller than v. The corresponding wavofunction repre- 
sents the infalling partner particles of the outgoing Hawking 
radiation, which have a negative energy as measured from in- 
finity. During the evolution, the high-wavenumber in- modes 
k"^ {x decreases ~» v increases) are being converted into the 
low-wavenumber Hawking radiation plus partner particles - 
where the breakdown of the JWKB approximation near the 
horizon leads to a mixing of these modes resulting in particle 
creation. 

For a stationary metric as in Eq. JQ), we may sepa- 
rate the most general solution of the Klein-Fock-Gordon 
equation |(2Jl into stationary modes with frequencies ui 
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= 0. 
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The black-hole horizon is assumed to be located at x = 
and hence the Taylor expansion of the velocity around 
this point reads 



(4) 



v{x) = -1 + KX + 0(K^a;^) 
with K denoting the surface gravity. 

B. Assumptions 



Let us summarize the assumptions that will be used 
for deriving the Hawking effect: 

Obviously, the surface gravity of the black hole (and 
hence the temperature of the Hawking radiation) must be 
much smaller than the cut-off scale where the concept of 
geometry and metric breaks down (i.e., F is not negligible 
anymore) 



K <m. fccutoff ■ 



(5) 



Furthermore, we shall assume that particle creation - ne- 
cessitating a break-down of the JWKB (geometric optics) 
approximation - occurs in the vicinity of the horizon only. 
Hence we shall consider an intermediate regime: close to 
the horizon at x = in units of n 

(6) 

but still many cut-off lengths away from the horizon 

I a; I fccutoff 1 ■ (7) 

Based on the above assumptions, we may neglect terms 
of second and higher order in k and a; since we are inter- 
ested in low-frequency modes uj = 0{k) only (Hawking 
radiation). Accordingly, the wave equation Q simplifies 
to 
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ox'^ ox 
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At this stage, it is advantageous to Laplace transform 
this equation via 



(9) 



where the contour C in the complex plane will be dis- 
cussed below. Note the change of sign in the second 
term due to the integration by parts. The wave equation 
for the Laplace transformed mode (puji^s) in terms of the 
complex variable s reads 



d 2 
2k — s — 2[iuj 
OS 



k)s (j), 



(10) 



In the following, we shall impose the following condi- 
tions on the dispersion relation Fls^]: 
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The dispersion relation F[s^] is assumed to be an 
analytic function of s^. 



• Hence it possesses a Laurent/Taylor expansion 

2n 

„ \ "-cutoff 



(11) 



where the (non-vanishing) coefficients a„ and the 
radius of convergence are supposed to be of order 
one - i.e., the dispersion relation does not depend 
on small quantities like w/fccutoff ■ 

• Furthermore, we assume a sub-luminal dispersion 
relation, i.e., in the rest frame, we have 



cLu 
dk 



< F[~fc2] < . 



(12) 



• Finally, we assume that asymptotically fc^ t oo, the 
dispersion relation is well separated from the line 
uj = k, i.e., the phase velocity does not approach 
unity 

F\-k'^] 

lim — < 1 lim \^ ^ ^Foo>Q. (13) 
fc^too k^ fc2too k^ 

Apart from these assumptions we do not need to specify 
the dispersion relation any further. For convenience, we 
shall choose units in which fccutoff = 1 and omit it in the 
following equations. 



III. ANALYTICAL DERIVATION 



As in the usual Fourier transform, we choose a contour 
that approaches infinity along the imaginary axis s = ik. 
In this case, the overall exponent in Eq. Q is purely 
imaginary and behaves for large \s\ — \k\ as 



exp < xs 



ds 



2k 



exp < xs 
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according to assumption H13|) . Hence the exponential 
function is rapidly oscillating at large \s\ and thus yields 
(again at large |s|) no contribution to the integral in 
Eq. jnj. The fc-integral over exp{ifc[a;— Foo/(2k)]} gives 
S{2kx — Foo) and hence vanishes since k\x\ <C 1. (As we 
shall see below, the same result can be obtained by de- 
forming the contour of integration in the complex plane.) 
From a physical point of view, this result is not very sur- 
prising since - given a non-trivial dispersion relation - 
one would not expect momenta which are much larger 
than the cut-off to contribute. (This expectation is how- 
ever false for F = 0.) The significant contributions will 
be found by the stationary phase method, or, after de- 
forming the contour, the saddle point method. 



B. Saddle- Point Method 

In order to apply the saddle point method, let us 
rewrite the Laplace transformation in Eq. Q as 



cl,^{x) = J dsg{s)e'=f^'l 
c 

with the two auxiliary functions 

5(s) = — : — , 



(17) 



(18) 



A. Complex Plane and Asymptotics 

After a separation of variables, the Laplace trans- 
formed wave equation l(Tn|) can be cast into the following 
form 

OS'^'H- 2ks- (14^ 

Up to an irrelevant pre-factor due to the integration con- 
stant, its solution reads 



exp I I ds 



2k 



(15) 



Since F[s^]/s^ is analytic, the Laplace transform (p^is) 
has a singularity at s = and a branch cut from s — 
to infinity - but no further singularities at finite values 
of s. We choose the negative real axis Q{s) = and 
3ff(s) < for the branch cut since this choice will be 
most convenient for deriving Hawking radiation ~ for an 
alternative choice, see Section Ivl 



and 



X \J 2k 

The saddle points s* of /(s) are determined by 
'df' 



ds 



(19) 



= ^ 2Ksix + F[si] = 0. (20) 



Many cut-off lengths away from the horizon \x\ ^ 1 
(but still k\x\ <C 1), we can approximate the integral in 
Eq. lO by the saddle-point expansion 



where a sum over multiple saddle points is implied with 
proper orientation. The next terms of the saddle-point 
expansion are suppressed by a factor of order 



5"(«*) ^0^1 



g{st:)xf"{s^) \xs^. 



(22) 
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and can be neglected if x is sufficiently large to overcome 
the smallness of s*, which depends on the small quantity 
Kx via Eq. 

Along the imaginary axis s» = ik^, the possible saddle 
points are given by 

F[-kl] ^2Kxkl^ (1 + 0(^2), (23) 

i.e., they exactly coincide with the solutions of the dis- 
persion relation (see Fig. ^ for A: ^ w 



kl-F[~kl] = {uj + vk,f^v^kl, 
since ui = 0{k) and ^ 1. 



(24) 




C. Contour for a; < 

Since the solutions of the dispersion relation and hence 
the saddle points depend on the sign of i.e., on which 
side of the horizon is considered, it is convenient to choose 
different contours in the complex plane for a; > and 
X < making sure that they are deformable to each 
other as x goes through zero. Let us first study the case 
a; < 0, i.e., the solution of the wave equation beyond the 
horizon, cf. Fig. [3 

Along the imaginary axis s — ik the exponent in 
Eq. H21(l is purely imaginary 



Im(s) 



n{f{tk)}^Q, 
whereas its derivative is purely real and positive 



rm = 1 - 



F[-k^] 
2k X fc2 



> 0. 



(25) 



(26) 



because x < and F[—k'^] > 0, cf. assumption (|12|l . 
As a result | exp{2;/(s)}| decreases rapidly {x < 0) with 
increasing 5R{s} and there are no saddle points on the 
imaginary axis s = ik. Hence we can deform the contour 
into the valley at 5R{s} > until we hit possible saddle 
points at 



with values 



^{s4x < 0]} > 0, 



^{f{s4x < 0])} > 0. 



(27) 



(28) 



Since the coefficients in the Laurent /Taylor expansion of 
F[s2] are of order one, the real part of the saddle point 
satisfying 2k.s1x + F[sl] = is mainly determined by 
the small quantity kx 1. Again we assume that the 
size of X ^ 1 overcomes the smallness of k a; ^ 1 and 
consequently obtain a solution 



< 0) 



27r 



s-«^-)g(s,), (29) 



which decays exponentially fast beyond the horizon. 



Re(s) 



FIG. 2: Landscape plot (top) of the real part of the logarithm 
of the integrand in Eq. @ for a sub-luminal dispersion rela- 
tion and the case x — —5 and ui = k — 1/30 as well as contour 
in the complex plane (bottom). The black square denotes the 
singularity and the zig-zag line is the branch cut. The behav- 
ior of the landscape near the imaginary axis is generic, but the 
structure away from that axis (e.g., existence of further saddle 
points) depends on the particular form of the (sub-Iuminal) 
dispersion relation (here F[s'^] = s*). 



Note that, even if we encounter no saddle points, the 
contour can still be deformed such that the contributions 
are exponentially small, cf. Fig. [3 In this way the cho- 
sen contour yields basically no contribution beyond the 
horizon x < - which is exactly what we want for the 
derivation of the outgoing Hawking radiation. 



D. Contour for a:: > 

In order to derive the solution outside the horizon x > 
0, another contour is needed for applying the saddle-point 
method, cf. Fig. 13 The exponent in Eq. (|21|l is still purely 
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FIG. 3: Landscape plot (top) of the real part of the logarithm 
of the integrand in Eq. @ for a sub-luminal dispersion rela- 
tion and the case x = +5 and u = k. — 1/30 as well as contour 
in the complex plane (bottom). The black square denotes the 
singularity, the black dots are the saddle points, and the zig- 
zag line is the branch cut. The behavior of the landscape 
near the imaginary axis is generic, but the structure away 
from that axis (e.g., existence of further saddle points) de- 
pends on the particular form of the (sub-luminal) dispersion 
relation (here F[s'^] ~ s'^). 



imaginary along the imaginary axis s 

F[-fc2] 



1 



ik, but the slope 



(30) 



2k X k"^ ' 

changes its sign at saddle points sf — ±ik^. The condi- 
tion ensures the existence of exactly two (symmetric) 
saddle points along the imaginary axis ~ i.e., solutions of 
the dispersion relation with finite values of /"(s*) - but 
the analysis can easily be generalized to the case of more 
than two saddle points. 

For small fc <C 1, the first term dominates according to 
Eq. (HH) 



F[-fc2 t 0] < /c^ 



F[-k^] 
2k X k^ 



(31) 



and hence the valley is on the side 5ft{s} < of the imag- 
inary axis - whereas for k^ > k^, the slope f'{ik) is neg- 
ative and thus the valley is on the other side ^{s} > 0. 
Hence the contour must cross these two saddle points 
(cf. Fig. O and pick up the corresponding contributions 



0±(x>O) 



2tt 



''fist) 



gist) 



2n 



-xf"{±ik^ 



2lT 



Txf'iik^) 



^±xf(ik,] 



g{±ik^:) 



gi±ik,). (32) 



(t)uist,x > 0) 




9(4) 




4>^{s^,x > 0) 




gis*) 





As f{ikf) is purely imaginary, the only difference in the 
absolute values of the two contributions is determined by 
the branch cut in g{s) 



(33) 



Ergo, the two saddle points at sf = ztik^, yield two 
rapidly {x 1) and oppositely oscillating contributions, 
whose absolute values satisfy the above relation (which 
will become important later on). 



E. Branch Cut 

Between the two saddle points —k^, < k < fc*, the 
valley lies on the same side 5R{s} < of the imaginary 
axis as the branch cut does, cf. Fig. 13 If it was not for 
the branch cut, the contour (originating from infinity) 
could be closed in this valley after crossing the two saddle 
points at = ±ifc* such that all additional contribu- 
tions (possibly further saddle points) are exponentially 
smaller than those of the saddle points at sf = ±ik^. 
However, the branch cut demands that we integrate along 
it to s = from both sides (with the proper orientation) 
and in circumventing the branch cut, we pick up the dif- 
ference in the values of g(s) 



/TTWN \S\ 



•J / K 



5(3(5)iO)-g(5(s)TO) = 2sinh 

\ K / S 

In this way, we obtain an additional contribution 
^Tix) = 2sinh(^) J^\sr-/^x 

2k s2 



(34) 



X exp < xs 



(35) 



where the integral runs from along the negative real 
axis up to the intersection point of the contour with the 
branch cut. In view of x 1, we can omit the second 
term in the integrand and extend the interval to — oo 



>0) 



2 sinh — 

V K 



ds exp {xs} 
s UK^/'' 



(36) 
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Outside the horizon a; > - the region we are interested 
in - we may substitute x— ^'^d obtain 

oo 

0r(x>O)«-2sinh(^).-/K| (37) 



which is just an integral representation of the F-function, 
i.e., 

C'*(a; > 0) » -2sinh (^) T (^-^) x'"'" ■ (38) 

Together with the contributions in Eq. H32|) . this com- 
pletes the (approximate) solution of the wave equation 
outside the horizon - for the case that the solution ba- 
sically vanishes beyond the horizon. Of course, we can 
only draw this conclusion if the two contours for x > 
and a; < in Figs. |2] and O are deformable to each other 
as X crosses zero. This property is ensured by assumption 
since the part of the complex plane covered during 
the deformation of the contours {k\x\ ^ 1 ^ |s*| <C 1) 
is well inside the radius of convergence of order one. 

F. Bogoliubov Coefficients 

Let us identify the various parts of the solution. The 
contribution generated by the branch cut (/)™''*(x) is 
the wavefunction of an outgoing particle with a low 
wavenumber k = 0{lo) (e.g., Hawking radiation). The 
saddle-point contributions (/)^(x), on the other hand, are 
rapidly oscillating, since the largeness of x ^ 1 is sup- 
posed to be stronger than the smallness of s,(kx). 

As one can observe in Fig. ^ the group velocity of 
the low-energy mode 4''^^^^{x) exceeds v{x), as one should 
expect for an outgoing particle - whereas the group ve- 
locity of the rapidly oscillating modes (t^tix) is smaller 
than v{x). Hence these are the in-modes 0™. 

Furthermore, the frequencies of the rapidly oscillating 
modes in the freely falling frame Vl± = w ± uA:* have dif- 
ferent signs because A;* a; (although w > for both 
modes). As a result, the low-energy outgoing particle 
(e.g.. Hawking radiation) is a mixture of positive and neg- 
ative frequency (with respect to the freely falling frame) 
in-modes - which can be described in terms of the Bo- 
goliubov coefficients 

= a^cf'^ + PJ'^ . (39) 

A non-vanishing Bogoliubov /3i^-coeSicient of course cor- 
responds to the phenomenon of particle creation - i.e., 
the in-vacuum with respect to the freely falling frame 
is converted into a quantum state containing particles 
(Hawking radiation) by the horizon. The ratio of the 
Bogoliubov coefficients is determined by Eq. 133|) 

1/3^1 =e-^'^/«|a„|, (40) 

which is the well-known relation leading to the thermal 
Hawking spectrum. E.g., applying the unitarity relation 



of the Bogoliubov coefficients a • a''' — /3 • = 1, we im- 
mediately obtain the thermal spectrum {Nui) = cx; 
1/(6^'^"/'^ - 1). 

Note that the quantum state generated by the time- 
evolution of the in-vacuum also contains particles with 
low wavenumbers beyond the horizon (see Fig. As 
required by energy conservation and unitarity, for each 
outgoing particle of the Hawking radiation, there is a 
corresponding partner particle with negative energy (as 
measured from infinity) inside of the black hole - but this 
does not alter presented calculation, cf. |3[. There are, 
however, correlations between the outgoing Hawking par- 
ticle and its partner beyond the horizon - which generate 
the true thermal character (mixed state instead of pure 
state) of the Hawking radiation for any outside observer 
(thermo-field formalism 4,^]), see Section Ivl below. 

The branch cut in the complex plane caused by the 
horizon turns out to be a main ingredient for deriving the 
Hawking effect - it generates the contribution cf)]^^^ {x) as 
well as the ratio in Eq. (|33|l - which are both essential 
features. 



IV. SUPER-LUMINAL DISPERSION 




FIG. 4: Super-Iuminal dispersion relation (not to scale). The 
points of intersection (black circles) with the two lines for 
V > 1 (i.e., a; < 0) and v < 1 (i.e., x > 0) determine the 
solutions of the dispersion relation for a given ui. As in the 
sub-luminal case, the out-modes are the low-wavenumber so- 
lutions corresponding to the outgoing particles of the Hawking 
radiation (x > 0) and their infalling partners (x <Q). How- 
ever, the high-wavenumber in-modes k"^ have group velocities 
exceeding v and hence are approaching the horizon from the 
inside a; < 0. 

So far, we restricted our attention to a sub-luminal dis- 
persion relation only. As we shall see now, the case of a 
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super-luminal dispersion as in Fig. 0] can be treated in 
basically the same way. The steps and derivations from 
Eq. Q to Eq. are identical, and we choose the same 
branch cut. Of com'se, for a super-luminal dispersion 
relation we have to modify assumptions (|12l) and H13() 
accordingly in order to ensure a vanishing asymptotical 
contribution at the imaginary axis s = ik. The deriva- 
tions in Sections IIII Al and IIII Bl apply in the same way, 
but now the solutions of the dispersion relation with large 
wavenumbers (the in-modes) are super-luminal, i.e., they 
originate from inside the black hole, see Fig. 0] 

A. Contour for a; > 

Let us first consider the solution outside the black hole, 
see Fig. In contrast to the sub-luminal case, the func- 
tion /(s) in the exponent has a positive slope for all de- 
values 

because F[—k'^] < 0, and, consequently, the valley is 
now situated at negative real parts of s. Deforming the 
contour into the valley, all contributions become expo- 
nentially small - but (again) we have to circumvent the 
branch cut. The contribution of the branch cut yields 
the same result as in Section IIII El Ergo, outside the 
black hole, we have only the outgoing Hawking particle 
- which is exactly what one would expect in the super- 
luminal case. 



B. Contour for x < 

For a; < 0, i.e., beyond the horizon, the slope f'{ik) 
changes its sign at the saddle points (i.e., the solutions 
of the dispersion relation with large wavenumbers) and 
the contour has to cross the imaginary axis picking up the 
saddle-point contributions, cf. Fig.lHI For those contribu- 
tions, we basically obtain the same results as in Eq. 1)32(1 
and hence in Eq. H33|) . because the branch cut is identical. 

Therefore, we reach the same conclusion as in Section 
IIII Fl where now the in-modes originate from inside the 
black hole. Hence we reproduce Hawking radiation also 
for a super-luminal dispersion - provided that the in- 
modes (large wavenumbers) are initially in their ground 
state with respect to the freely falling frame. 

V. ENTANGLEMENT 

So far, we restricted our attention to the decomposi- 
tion of the outgoing Hawking radiation in terms of the 
in-modes, cf. Eq. (|39|l . However, a full description of the 
evolution of the quantum state requires a complete set of 
out- modes, i.e., the outgoing Hawking particles as well as 
their infalling partners. Fortunately, it turns out that the 




FIG. 5: Landscape plot (top) of the real part of the logarithm 
of the integrand in Eq. JUJ for a super-luminal dispersion rela- 
tion and the case x = +5 and ui = k, = 1/30 as well as contour 
in the complex plane (bottom). The black square denotes the 
singularity and the zig-zag line is the branch cut. The behav- 
ior of the landscape near the imaginary axis is generic, but the 
structure away from that axis (e.g., existence of further saddle 
points) depends on the particular form of the (super-luminal) 
dispersion relation (here F[s'^] = — s*). 



Bogoliubov coefRcients of the infalling partners can be in- 
ferred in complete analogy to the previous Sections if we 
choose the branch cut in the opposite way, i.e., along the 
positive real axis 3(s) = and ^{s) > 0. Circumvent- 
ing this branch cut then reproduces the wavefunction of 
the infalling partner particles as in Sec. IIII El (but with 
X — > —x), and, consistently, this contribution only occurs 
for a; < 0, i.e., beyond the horizon. Again, for both cases 
(sub- and super-luminal) , one obtains basically the same 
relation 

^partner ^ ^partner ^in _^ ^partner ^in ^ (-42) 



8 
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Re(s) 



FIG. 6: Landscape plot (top) of the real part of the logarithm 
of the integrand in Eq. JUJ for a super-luminal dispersion re- 
lation and the case x — —5 and ui = k = 1/30 as well as 
contour in the complex plane (bottom). The black square 
denotes the singularity, the black dots are the saddle points, 
and the zig-zag line is the branch cut. The behavior of the 
landscape near the imaginary axis is generic, but the struc- 
ture away from that axis (e.g., existence of further saddle 
points) depends on the particular form of the (super-luminal) 



dispersion relation (here F[s 



Note that we have interchanged the role of the cre- 
ation and annihilation operators and hence Q-P^-ftner 
^partner j^gj.^ bccause the energy of the infalling partner 
particles is negative as measured from infinity and their 
pseudo-norm is negative for positive lo. The opposite 
direction of the branch cut implies the inverse relation 
compared to Eq. H33|l and together with the above inter- 
change, we derive the same ratio as in Eq. (|40|l 



1/3, 



partncri —7ra;/K i partner i 



(43) 



The knowledge of the complete set of out-modes (Hawk- 
ing radiation ^Hawkmg pj^g their infalling partners 



^partner^ facilitates the decomposition of the in-modes in 
terms of the out-modes 



(44) 



In view of the relations (|40|l and H43|l as well as unitarity 
a ■ a"^ — (3 ■ (3'' — 1, the inverse Bogoliubov coefficients 
satisfy an analogous condition 



1/3, 



invi _ g-Troj/Kl^invl 



(45) 



Consequently, the in- vacuum defined via ain |0i„) = will 
be annihilated by a linear combination of the operators 
corresponding to the out-modes 



'-partner _|_ — ttcij/k / - Hawking\ f 



|0i„ 



0, 



(46) 



where an irrelevant phase has been absorbed by the re- 
definition of dP'*''*"°''. This well-known relation (see, e.g., 
^0]) induces the entanglement between the particles of 
the Hawking radiation (jHawking their infalling part- 
ners aPartncr 

|0)r„ (X exp (aP-tncr -Hawking) 1 1 j^^^^^ ^^^^ 

which in turn generates the thermal density matrix after 
averaging over the unobservable infalling partners. 



VI. UNIVERSALITY 

The derivation presented in the previous Sections 
demonstrates that the the Hawking effect does (to low- 
est order) not depend on the details of the dispersion 
relation at high wavenunibers - given the model assump- 
tions discussed above. Let us try to identify more general 
conditions under which the Hawking should remain un- 
changed by the details of the physics at large wavenuni- 
bers. For convenience, we shall assume that the cut-off 
scale coincides with the Planck scale and use the terms 
sub-Planckian for effects according to the known laws of 
physics (e.g., linear dispersion) and trans-Planckian for 
new physics (e.g., non- linear dispersion). 

First of all, we assume that the JWKB (geometric op- 
tics or eikonal) approximation breaks down (thereby al- 
lowing for the phenomenon of particle creation) in the 
vicinity of the horizon only, where the gravitational red- 
shift induces a transition of trans-Planckian into sub- 
Planckian modes. An example where this assumption 
does not apply will be discussed in Section IVl Bl 

Given that assumption, the crucial point is the quan- 
tum state of the modes when they leave the Planckian 
regime. If the modes leave the Planckian regime ("are 
born" ) in their ground state with respect to freely falling 
observers near the horizon, then one obtains Hawking 
radiation, cf. 

MM Let us review the standard ar- 
gument leading to that conclusion. In terms of the 
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Regge- Wheeler tortoise coordinate r*, the 1+1 dimen- 
sional Schwarzschild metric can be cast into the confor- 
mally flat form 



(48) 



where the ~ applies near the horizon. The trajectory of a 
freely falling observer is given by [A and B are integration 
constants) 



r,(i t oo) ~ - ^ exp <j \+B. 



2MJ 



(49) 



and its proper time dr^ = ds'^[t,r^{t)] accordingly reads 



T ~ exp 



2M ] 



(50) 



Hence the freely falling observers would define their 
ground state via the positive frequency solutions 



1 



^ — iioU 



(51) 



with respect to the Kruskal coordinate 

U = -AM e-"/(4^^) = -AM e"[*-'-*l/(4*^) . (52) 

The doubly exponential behavior of these modes - when 
expressed in terms of the coordinates t, r of an outside 
observer - lead to the thermal particle content. 

The remaining issue is, of course, to determine in which 
cases the modes do indeed leave the Planckian regime in 
their ground state with respect to freely falling observers 
(near the horizon). As a very natural example, one could 
ensure this property by means of the following three as- 
sumptions: 

a) Freely falling frame 

If we assume that the usual local Lorentz invariance 
is broken at the Planck scale via the introduction 
of preferred frames (where preferred frames are the 
frames in which Planckian physics displays maxi- 
mal symmetry under time-inversion, for example) 
then the freely falling frame should be preferred 
(instead of the rest frame of the black hole, for ex- 
ample) . 

b) Ground state 

The Planckian excitations are assumed to start off 
in their ground state (with respect to the freely 
falling frame, see point above) subject to possible 
constraints such as conservation laws etc. 

c) Adiabatic evolution 

Finally the evolution of the modes is supposed to 
be adiabatic - i.e., the Planckian dynamics is sup- 
posed to be much faster than all external (sub- 
Planckian) variations (e.g., experienced by a trav- 
eling wavepacket). This condition demands the 



absence of level crossing and long time-scales in 
Planckian physics. 

E.g., for the sonic black hole analogues ("dumb holes") 
such as a fluid flowing trough a Laval nozzle (accelerated 
from sub-sonic to super-sonic speed), the freely falling 
frame corresponds to the local rest frame of the flow- 
ing fluid, whereas the rest frame of the walls of the noz- 
zle is analogous to the global rest frame of black hole. 
Of course, one can easily imagine situations where at 
least one of the above assumptions fails. E.g., for a 
super-luminal dispersion relation, the modes with large 
wavenumbers originate from inside the black hole, i.e., 
ultimately from the singularity (of from a turbulent 
regime), and it is not obvious why they should be in 
their ground state. Further examples for the failure of the 
above assumptions, where the reference frame for Planck- 
ian physics is not the local freely falling frame but the 
global rest frame of black hole; or where the adabaticity 
breaks down, are the subject of the next Sections. 



A. Miles Instability 

As an example, in which the aforementioned set of as- 
sumptions fails and which does not reproduce Hawking 
radiation, let us consider the following fluid model: Apart 
from a deviation from the usual dispersion relation as de- 
scribed by a /c-dependent phase velocity Wph(fc); we sup- 
pose a coupling to a reservoir of Planckian degrees of 
freedom in the rest frame of the black hole (i.e., not the 
freely falling frame) manifesting itself as an effective dis- 
sipation term in the dispersion relation 



{U! + Wflfc)^ = k 



(53) 



The fc-dependence of the damping term 7(fc) ensures that 
it is completely negligible at sub-Planckian wavenumbers. 

In flat space-time (fluid at rest vr = 0), the damping 
term just implies a decay of the Planckian modes (as one 
would expect). For a black hole, the Planckian modes 
giving rise to Hawking radiation, however, behave in a 
different way. As one can easily perceive from Figs. ^ 
and 01 for solutions of the dispersion relation with large 
(Planckian) wavenumbers, we have 



hph(^) 



< 1 



(54) 



Expanding the relevant solution of Eq. H53|) for uj in pow- 
ers of this small quantity, we obtain 



2 uflfc + i^{k) 



(55) 



Assuming a real wavenumber, the imaginary part of the 
frequency changes its sign if the fluid velocity exceeds the 
phase velocity 



3(w) > 0, 



(56) 



10 



which indicates an instabiHty. This phenomenon is basi- 
cally the Miles instability - which is responsible for the 
generation of water waves by wind, for example |l2j . 

If k and 7 are of order one (in Planckian units), the 
imaginary part of lu is of the same order as the real part 



(57) 



and since uj = 0{k) corresponds to the inverse size of the 
black hole, there can be enough time for the instability 
to develop and to excite the modes. Note that positive 
frequency (trans-Planckian) modes with Vp^{k) > v"^ are 
damped but negative frequency modes with v^Yiik) < w| 
are amplified. Hence this effect destroys the balance in 
Eq. (|40|l which generates the thermal spectrum of the 
Hawking radiation. 

However, the above analysis based on classical solu- 
tions of the dispersion relation cannot be applied directly, 
i.e., without respecting the fluctuation-dissipation theo- 
rem, for example, to the quantum fluctuations that gen- 
erate the Hawking radiation. In order to turn our atten- 
tion to the quantum theory, let us consider the following 
Lagrangian density corresponding to a super-luminal dis- 
persion Wpjj = 1 + 



(V^ 



(V^ 



(58) 



For a stationary metric, i.e., v{x), a conserved energy 
density with respect to global rest frame of black hole 
can be derived by means of the Noether theorem 



(V<i 



(59) 



Evidently the energy density is not positive definite for 
> Wpj^, i.e., beyond the horizon (super-luminal dis- 
persion). The local energy density with respect to the 
local freely falling frame is of course positive definite. 
After a normal mode expansion into wavepackets, the 
total Hamiltonian is split up into nearly independent pos- 
itive and negative energy modes (with respect to global 
rest frame of black hole). Obviously, the negative energy 
modes can be strongly excited by a comparably weak in- 
teraction with further Planckian degrees of freedom at 
the global rest frame of the black hole. In this way, these 
modes would not be in their ground state - even with re- 
spect to freely falling observers - and, consequently, one 
would not reproduce Hawking radiation. 

For example, let us consider fermionic fields where the 
quantum states of all trans-Planckian modes are maxi- 
mally excited, i.e., |1) instead of the ground state |0). In 
that situation, the usual relation /J^^ = e"'^"/''^^; implies 



(iliV^ID - 



2 _j_ ^ — 2tt(jj/k. 



(60) 



i.e., not a thermal spectrum. Note that, for fermions, 
the unitarity relation is a • _|_ ^ . _ instead of 
a • a''' — /3 • = 1 leading to the Fermi-Dirac spectrum 



-[^^^g27ra)/K _^ ^Yie usual Hawking radiation. If the 

occupation number (i.e., |0) or |1)) depends on the his- 
tory of the mode - e.g., the frequency lu - then one would 
obtain another (in general non-thermal) spectrum. With 
an appropriate mixture of |0) and |1), one could even ob- 
tain a state with no outgoing particle content (Hawking 
radiation) at all (Boulware vacuum). 

In summary, an interaction with a reservoir at the 
Planck scale with respect to the rest frame of the black 
hole can invalidate Hawking's derivation. In that argu- 
ment, the rest frame of the black hole is a crucial point - 
a damping term with respect to the freely falling frame 
2i(uj + Vfik)j(k) would not induce a positive imaginary 
part of w. A similar phenomenon occurs in the so-called 
" black hole laser" where wavepackets bounce back and 
forth between the inner and outer horizons - which also 
generates deviations from the Hawking effect 13:|. This 
quasi-reflcction mechanism also singles out the rest frame 
of the black hole (location of the two horizons) as a pre- 
ferred frame for the Planckian modes. In contrast to the 
Miles instability, this phenomenon displays more similar- 
ities to the Pierce instability . 



B. Breakdown of Adiabaticity 



In the previous subsection IVI Al the assumption a) of 
Section IVll and hence also b) failed. Let us now give an 
example for the breakdown of the adiabaticity condition 
c), which is closely related to the assumption that geo- 
metric optics is valid everywhere except in the vicinity of 
the horizon. 

One version of the adiabatic theorem states that if the 
dynamics of all internal degrees of freedom is much faster 
than any external time dependence, then a system being 
initially in its ground state basically remains in the (time- 
dependent instantaneous) ground state. (Of course, for 
this theorem to apply we have to assume that quantum 
theory is still valid at the Planck scale.) 

As a counter-example, where the system does not stay 
in its ground state, consider the dispersion relation 



sin^ k + 



(61) 



which has minima at fc € ttN (Planckian units). If we 
assume a weakly time-dependent metric far away from 
the black hole ds^ = a^{t)[dt'^ — dx^], the wave equation 
reads after a normal mode expansion 



^fc -I- [sin^ k + a^{t) m^] <j)k ^ , 



(62) 



since the mass term breaks the conformal invariance. 
Therefore, it very easy to create Planckian (fc sa ttN) 
particles (e.g., via parametric resonance) by means of 
comparably small and slow (sub-Planckian) variations of 
a{t) with with a characteristic scale corresponding to m 
(instead of the Planck mass). 

As a result of this breakdown of the geometric optics 
approximation far away from the horizon, the Planck- 
ian modes falling towards the black hole are not in their 
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ground state - and hence one will again obtain deviations 
from the Hawking effect. Note that a similar effect (oc- 
cupation of Planckian modes) can occur during inflation 
if we assume a dispersion relation like the above. 

VII. CONCLUSIONS 
A. Summary 

The Hawking effect is not a priori independent of the 
laws of physics at the Planck scale, but it can be made 
so by imposing the three assumptions a) Freely falling 
frame, b) Ground state, and c) Adiabatic evolu- 
tion, explained in more detail in Section IVll As one ex- 
ample, we generalized the analytical method of Ref. Q 
to arbitrary dispersion relations subject to some rather 
general assumptions. 

However, we have also demonstrated counter- 
examples, which do not appear to be unphysical or artifi- 
cial, displaying deviations from Hawking's result. There- 
fore, whether real black holes emit Hawking radiation or 
not remains an open question and gives non-trivial infor- 
mation about Planckian physics. 

B. Outlook 

Another example, where sub-Planckian phenomena 
have their origin in trans-Planckian modes, is the gen- 



eration of inhomogeneities during the cosmic epoch of 
inflation (according to our present standard model of cos- 
mology) from quantum fluctuations of the inflaton field. 
In this case, the investigation of the universality, or, con- 
versely, the dependence of this mechanism on Planckian 
physics including higher-order corrections has an addi- 
tional aspect, because observations of the cosmic micro- 
wave background, for example, might yield signatures of 
Planckian physics, see, e.g., [laillfl- E.g., a dispersion 
relation with minima for large fc- values as in Section fVl HI 
potentially allows par ticle creation leading to a change in 
the spectrum, see jT^- 
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